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Abstract 



' O We employ a coarse-graining approach to analyze nonlinear cascades in Boussinesq flows using 

high-resolution simulation data. We derive budgets which resolve the evolution of energy and poten- 
tial enstrophy simultaneously in space and in scale. We then use numerical simulations of Boussinesq 
flows, with forcing in the large-scales, and fixed rotation and stable stratification along the vertical 
Y^ axis, to study the inter-scale flux of energy and potential enstrophy in three different regimes of 

■ ,—( stratification and rotation: (i) strong rotation and moderate stratification, (ii) moderate rotation 

^ and strong stratification, and (iii) equally strong stratification and rotation. In all three cases, we 

observe constant fluxes of both global invariants, the mean energy and mean potential enstrophy, 
from large to small scales. The existence of constant potential enstrophy flux ranges provides the 
first direct empirical evidence in support of the notion of a cascade of potential enstrophy. The 
persistent forward cascade of the two invariants reflects a marked departure of these flows from 
Cn two-dimensional turbulence. 
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1 Introduction 

The Kolmogorov [T] and Kraichnan [2 theories of three- and two-dimensional Navier-Stokes turbulence 
have served as a benchmark in the understanding of fluid turbulence and as fundamental tests for the 
accuracy of simulations. The Boussinesq approximation of the compressible Navier-Stokes equations in 
a rotating frame give a fairly accurate description of the flow dynamics over much of the Earth's oceans 
and atmosphere but are prohibitively expensive to simulate in detail over global scales. Guided by the 
success of the Kolmogorov/Kraichnan theories, it would be useful to develop a statistical phenomenology 
of the small scales of Boussinesq flows to gain an understanding of the physics, serve as benchmarks to 
test simulations, and offer parameterizations which could eventually be useful in practical modeling of 
geophysical flows. In addition to global energy, the inviscid Boussinesq equations conserve local potential 
vorticity (hereafter, PV) and global potential enstrophy, thus offering more complexity than incompress- 
ible Navier-Stokes dynamics. Charney [3 addressed one limiting case of the Boussinesq approximation, 
namely the quasi-geostrophic limit of strong rotation and strong stratification, and showed that the 
conservation of both energy and the quadratic potential enstrophy in such flows constrained energy to 
cascade to the large scales as in 2D turbulence [Ij. 

Conservation of potential enstrophy is believed to play a fundamental role in the dynamics of the 
atmosphere and oceans [4]. Understanding its function in nonlinear scale interactions would appear 
to be essential for extending Kolmogorov's theory to Boussinesq flows with rotation and stratification. 
Herring et al. [5j studied the cascade properties of potential enstrophy in turbulence simulations with a 
passive scalar in the absence of rotation and stratification. The authors concluded that because potential 
enstrophy in their simulation is not a quadratic, but a quartic invariant, the usual Kolmogorov-like 
arguments for a cascade and an inertial range do not apply. In fact, their potential enstrophy budget 
shows the direct action of viscous-diffusion terms at all scales, even in the limit of very small viscosity 
and diffusivity. Their figure 15 shows that potential enstrophy dissipation peaks at the largest scales 
(see also their figure 14 and the discussion on pp. 37 and 43). Thus, a pure inertial range of potential 
enstrophy flux is precluded because of contamination by dissipation at all scales. However for Boussinesq 
flows, with strong rotation and/or strong stratification, Kurien et al. jB] (hereafter, KSW06) derived 
analytically, starting from the evolution of the two-point correlation of potential vorticity, a flux law for 
potential enstrophy which is analogous to Kolmogorov's 4/5-law for energy flux in 3D incompressible 
turbulence. The so-called 2/3-law of KSW06 implies that an inertial cascade of potential enstrophy can 
exist in three limiting cases: (i) strong rotation with moderate stratification (ii) moderate rotation with 
strong stratification and (iii) strong rotation and strong stratification. By "moderate" we mean that the 
rotation (stratification) frequency is of the same order as the non-linear turbulence frequency; whereas 
by "strong" we mean that the respective frequencies are much faster than the non-linear timescale. In 



the above three regimes, KSW06 showed that potential cnstrophy, generally a quartic quantity, becomes 
quadratic which results in the localization of viscous-diffusion terms to the smallest scales, thus allowing 
for an inertial range of scales dominated by the flux of potential enstrophy. Furthermore, KSW06 
suggested that in the absence of strong rotation and/or strong stratification, when potential enstrophy 
reverts to being quartic, viscous-diffusion effects may contaminate all scales, in agreement with the 
conclusions of [5]. 

The existence of an inertial cascade range for potential enstrophy is far from obvious and remains an 
unsettled issue. To date, there has been no empirical demonstration of KSW06's results on the constant 
flux range of potential enstrophy in Boussinesq flows. In a phenomenology and supporting data from 
Boussinesq simulations with equally strong (non-dimensional) rotation and stratification were presented 
to show that conservation of quadratic potential enstrophy could constrain the spectral distribution 
of energy in the large wavenumbers. In [5] the analysis and phenomenology of 7] was extended to 
include the two other limiting cases (i) and (ii) above, to show that indeed in all parameter regimes with 
linear PV, the downscale flux of (quadratic) potential enstrophy can constrain the scale-distribution of 
energy. While the two studies [71 [S] did not measure potential enstrophy flux, their numerical findings 
confirmed phenomenologically predicted scalings of the energy spectra assuming constant downscale 
fluxes of potential enstrophy. Motivated partially by those results, the present work uses data from [8] 
to directly measure the cascade of potential enstrophy for the first time. 

Some previous studies have highlighted the interest in energy fluxes. The results from HO] of purely 
stratified flow computed in small-aspect-ratio domains with forcing of the horizontal velocity only, did 
not show convincing constant energy fluxes (see figure 18 of [5] and figure 15 of [TD] which plot energy 
flux on a log- log scale). On the other hand |TT] provides evidence of scale- independent fluxes of kinetic 
and potential energy in purely stratified flow in unit aspect-ratio. It should be noted that these studies 
used forcing, stratification, and domain aspect-ratio different from ours. Furthermore, the simulations 
of [HI IIOI 111] were not reported to be in the linear PV regime and, therefore, our results may not apply 
to their flows. Most importantly, none of these previous studies computed or analyzed the potential 
enstrophy flux, which constitutes an essential part of our work. 

In this Letter, we present a very general framework for analyzing nonlinear scale interactions in 
Boussinesq flows. The coarse-graining approach we utilize allows for probing the dynamics simultane- 
ously in space and in scale. Motivated by the work of [3 IH]j we then measure fluxes of energy and 
potential enstrophy across scales from simulations in three distinct limits of rotation and stratification. 
Our results show constant and positive fluxes of the two quadratic invariants, indicating simultaneous 
persistent downscale cascades of both quantities in all three cases. Our measurements of potential en- 
strophy flux are a novel contribution of this Letter and constitute the first empirical confirmation of 



analytical results by "5^. Furthermore, our evidence of a scale-independent energy flux is significant be- 
cause it conveys that a cascade should persist to arbitrarily small scales at asymptotically high simulation 
resolutions. 

2 Boussinesq dynamics 

We study stably stratified Boussinesq flows in a rotating frame. The dynamics is described by momentum 
^ and active scalar ([2| equations: 

dtu + (u-V)u = -Vp - /zxu - N9z + j/V^u + :F", (f) 

dt0 + iu-V)e = Nu^ + K,V^e + J-'^. (2) 

Here, u is a solenoidal velocity field, V-u = 0, whose vertical component is Uz. The effective pressure is 
p, and ^F^ ,J^^ are external forces. Gravity, g, is constant and in the — z direction. Total density is given 
by Pt(x) ^ Po — bz + p(x), such that |p(x)| ^ \bz\ and |p(x)| <C Pq, where po is a constant background 
density, b is constant and positive for stable stratification, and p(x) is the fluctuating density field with 



zero mean. The normalized density, 0(x) — y^g/bpop{x), has units of velocity. For a constant rotation 
rate 51 about the z-axis, the Coriolis parameter is / = 2il. The Brunt- Vaisala frequency is A^ = y'gb/pQ, 
kinematic viscosity is i/ = m/POj and mass diffusivity is k. In this paper, we only study flows with Prandtl 
number Pr — v / k, — 0(1). Relevant non-dimensional parameters are Rossby number, Ro — fni/f, and 
Froude number, Fr = fni/N, where we define the characteristic non-linear frequency as fni = (efk'jy'^, 
for a given energy injection rate e/ at wavenumber kf (see |12l [5ll5]). 

The dynamics of inviscid and unforced Boussinesq flows (such that v = k = J'^ = Ff = 0) is con- 
strained by the conservation of potential vorticity, q{x) ~ -j-Wq-Vpt, following material flow particles, 
Dtq — dtq + {u-'V)q = 0. Here, absolute vorticity is Wq = a; + /z and local vorticity is a; = Vxu. PV 
may be written in terms of u; and 6 as 

q{x) ^ fd^e-Nuj-i + u^-Ve- fN. (3) 

The first two terms are linear and dominate over the quadratic term, u^-Vfl, in the limit of large / 
and/or large N. The constant part in (|3| does not participate in the dynamics and can, therefore, be 
neglected jl2j. 

In addition to conservation of PV, the flow is constrained by the global conservation of potential 
enstrophy, Q — \q^ , such that ^(Q) — 0, where (...) — -^ /^ d'^x(. . . ) is a space average. Another 
quadratic invariant of the inviscid dynamics is total mean energy, Et = |(|up + |6'p), such that ^Et = 0. 



3 Numerical data 

The Sandia-LANL DNS code was used to perform pseudo-spectral calculations of the Boussinesq equa- 
tions ([l|-(l2| on grids of 640'^ points in unit aspect-ratio domains. The time-stepping is 4th-order Runge- 
Kutta and the fastest linear wave frequencies are resolved with at least five timesteps per wave period. 
The diffusion of both momentum and density (scalar) is modeled by hyperviscosity of laplacian to the 
8th-power. The coefficient of the hyperviscous diffusion term is chosen dynamically such that the energy 
in the largest wavenumber shell (smallest resolved scale) is dissipated at each time-step [T31 [H] . This 
choice ensures that one does not have to guess the coefficient a priori and is a way to allow the flow itself 
to determine the magnitude of the diffusion. Hyperviscosity is a standard dissipation model that has 
long been used in studies of rotating and/or stratified flows O [101 HH IE] ■ In principle, hyperviscosity 
can lead to thermalization and isotropy at the smallest scales [^, however our results on the inertial- 
range cascades are robust and unaffected by the small-scale dissipation model ^TS\ . Stochastic forcing is 
incompressible and equipartitioned between the three velocity components and 6. The forcing spectrum 
is peaked at fcy = 4 ± 1, for large scale forcing. We use the two-thirds dealiasing rule. These data were 
reported in [71 [5], where further computational details may be found. 



Run resolution / (Ro) N (Fr) 



Rs 64(F 3000 (.002) 14 (.4) 

rS 640^ 14 (.4) 3000 (.002) 

RS 640^ 3000 (.002) 3000 (.002) 



Table 1: Parameters of the Boussinesq simulation data. 

We analyze three sets of simulations corresponding to three extreme flow regimes summarized in 
Table [I] The flrst, Rs, is a flow under strong rotation and moderate stratiflcation, f/N ^ 1. The 
second, rS, is a flow under moderate rotation but strong stratification, f/N <C 1. The third, RS, is a 
flow under strong rotation and strong stratification such that f — N. Figure [l] shows that in all three 
cases, (Q) is well approximated by (one half) the square of the corresponding linear PV to within 3% 
or better (see P, '5^ ) . Figure fl] is evidence that our simulations are in regimes of strong rotation and/or 
strong stratiflcation. We analyze snapshots of the flow at late times when (Q) along with small-scale 
energy spectra (at wavenumbers A: > 6) have reached a statistically steady state. The total energy, 
however, continues to grow due to an accumulation at the largest scales. 

4 Analyzing the cascades by coarse-graining 

Following [ini |T71 [ini m] , we use a simple flltering technique common in the Large Eddy Simulation 
(LES) literature to resolve turbulent flelds simultaneously in scale and in space. Other decompositions. 
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Figure 1: Time-series of mean potential enstrophy (Q) (solid line), and its quadratic part (Q) (crosses). 
In run RS (Q) = ( [fd^0 - Nco-zf )/2; in run rS (Q) = ( [Nu:-zf )/2; and in run Rs (Q) = ( [fd^Of )/2. 
The plots show that (Q) reaches steady-state and that (Q) ~ (Q) in all three regimes considered. 



such as wavelet analysis, also allow for the simultaneous space-scale resolution and may be used to 
analyze nonlinear scale interactions as well. 

We define a coarse-grained or (low-pass) filtered field in d-dimensions as 



a,(x)- /"d'*rG,(r)a(x + r), 



(4) 



where G(r) is a normalized convolution kernel, J d'^r G(r) = 1. An example of such a kernel is the 
Gaussian function, G{r) = -hpeT"^ ^^- Its dilation G^(r) = i^'^Givjt) in d-dimcnsions has its main 
support in a ball of radius I. Operation (H) may be interpreted as a local space average. In the rest of 
our Letter, we shall omit subscript i whenever there is no ambiguity. 

Applying the filtering operation Q to the dynamics ([l])-(|2]) yields coarse-grained equations that 
describe the evolution of U£(x) and ^^(x) at every point x in space and at any instant of time: 



^tu + (u-V)u 
d^ -K (u-V)5 



-Vp- fzxu- Ndz 

Nu,-V-T{u,e) 



(5) 
(6) 



where subgrid stresses, Ti{f,g) — fg^ — f^ g^, are "generalized 2nd-order moments" [17j accounting for 
the infiuence of eliminated fluctuations at scales < £. 



The coarse-grained equations describe flow at scales > ^, for arbitrary L The approach, therefore, 
allows for the simultaneous resolution of dynamics both in scale and in space. Furthermore, the approach 
admits intuitive physical interpretation of various terms in the coarse-grained balance (pi)- ([6]) which 
resemble the original governing eqs. (Il])-([2|, except for additional subgrid terms which quantify the 
nonlinear coupling between resolved and filtered scales. These subgrid terms depend inherently on the 
unresolved dynamics which has been filtered out. Traditional modeling efforts, such as in LES (see 
for example |20)). focus on devising closures for such terms which are plausible but whose regimes of 
applicability and validity are inevitably unknown. A key feature of the formalism employed here that 
distinguishes it from those modeling efforts is that it allows us to estimate the contribution of subgrid 
terms as a function of the resolution scale (, through exact mathematical analysis and direct numerical 
simulations (see for example [HJ [551 1131 IM])- Our approach thus quantifies the coupling that exists 
between different scales and may be used to extract certain scale-invariant features in the dynamics. 

4.1 Large-scale energy budget 

From eqs. ^ and ([6]), it is straightforward to derive an energy budget for the large-scales, which reads 



Here, i^jVup -I- k|V6'P is molecular dissipation acting on scales > £, and u-jF^ +6 F^ is energy injected 
due to external stirring. The term J£(x) represents space transport of large-scale energy whose complete 
expression is deferred to an Appendix below. Subgrid scale (SGS) flux, H^ (x), accounts for the nonlinear 
transfer of energy from scales > i? to smaller scales: 

n^(x) = -djUif{ui, Uj) - dj9 t{6, uj). (8) 

The SGS flux in ([8| is work done by large-scale velocity and scalar gradients, Vu(x) and V6'(x), against 
subgrid stresses. It acts as a sink in the large-scale budget ^ and accounts for the energy transferred 
across scale £ at any point x in the flow. Furthermore, Tlei^) is Galilean invariant. Other definitions of 
a flux are possible, such as Tli{x) — UiUjdj{ui — Uj) (Eq. (2.52) in Frisch [25j ). which differs from our 
definition (Is]) by a total gradient (disregarding the scalar part). However, these alternate definitions are 
not pointwise Galilean invariant, so the amount of energy cascading at any point x in the fluid according 
to such definitions would depend on the observer's velocity. 

Another physical requirement on the SGS flux H^ (x) is that it should vanish in the absence of fluc- 
tuations at scales smaller than t |211 122| . For example, when i = K^\^, where K^ax is the maximum 
wavenumber in a pseudospectral simulation, there should be no cascade across (, simply because fluctu- 
ations at wavevectors |k| > Kmax have zero amplitude. This is satisfied by our definition ([8]) identically 



at every point x in the flow. Alternate definitions such as nf(x) = UiUjdjUi (eq. (6) in '^S') fail this 
pointwise requirement of a flux. 

4.2 Large-scale Q budget 

Similar to the momentum, scalar, and energy equations ([5|-([7]), we can write down large-scale balances 
for PV and Q. The "bare" PV equation with diffusion and external forcing may be derived from eqs. 
and ^ as 

dtq + (u- V)g = v{Ve - Nz) -W^uja + nuia-^^Ve (9) 

+ {Ve ~ Nz) -T"^ + u;<,-VJ'^ 

where JF" = VXJF". Applying the filtering operation to ^k in the limit of strong rotation and/or 
stratification (limit of large / and/or iV such that PV is linear, q = fdzO — Nuj^), yields the following 
balance for large-scale PV 

dtq + (u- V)q = - V-r(g, u) - i^ N V^cJ^ + Kf V^dJ 

-NT^, + fd,T^. (10) 



Using (10), we can now write the large-scale potential enstrophy budget in the limit of linear PV, 



d,\sl + V.3Q = -UQ-DQ + eZ,, (11) 

where J'^(x) is space transport, D^{x) is dissipation due to viscosity and diffusivity acting directly on 
scales > £, and e^^ is the potential enstrophy injected due to external forcing. These terms are defined 
in [15] . n^ (x) is SGS fiux of potential enstrophy, 

Ilf{x) = -d,qr{q,u,). (12) 

It is straightforward to verify that ny(x) is Galilean invariant and vanishes in the absence of subgrid 
fiuctuations. 

5 Calculating fluxes using sharp spectral fllter 

We choose the so-called "sharp spectral filter" as our coarse-graining kernel in the definition of fluxes. We 
denote a field in a periodic domain [0, 1)"^, coarse-grained with the spherically symmetric sharp-spectral 
filter to retain only Fourier modes |k| < K, by 

a<-^(x)= J2 dka{k)e'^''^''. (13) 

|k|<_ft: 



This is similar to a^(x) with £ ^ K ^. We omit the factor 2tt in reference to wavenmnber in this Letter. 



While an isotropic filter such as ( 13 ) cannot distinguish between different directions, it does not average 
out anisotropy if present in a flow as will be clearly demonstrated in |15j . 

Using this filter, we can discern the amount of energy and potential enstrophy cascading across a 
certain wavenumber K. For example, to analyze the energy cascade, we can compute 

n;^(x) - -a,u<^T<^(«„«,) - 9jr^T<^(0,«,) (14) 

as a function of K. Here, T^^{f,g) — (fg)'^^ — f^^g^^ ■ We can also analyze potential enstrophy 
cascade using 

Ul{^)^-d,q<''T<^{q,u,). (15) 



SGS energy flux (14) coincides with that used in [21 [10] only after space averaging. Yet, our quantity 
has the correct pointwise physical properties discussed above and, therefore, allows for studying spatial 
properties of the cascades. 

6 Numerical Results 

In this Letter, we restrict our numerical investigation to spatially averaged fluxes using the sharp spectral 
filter. Figure [2] shows that there is a positive and constant flux (y-axis shown on a linear scale to highlight 
true constancy) of total energy to small scales in all three cases of rotation and stratificatiorR While 
all three cases show a clear downscale energy cascade, the RS case also has a negative flux over K < 4 
possibly indicative of the expected inverse cascade present in such regime (e.g. [271 [5S1|53] ) . We do not 
have the required scale-range in our simulations to say anything more definitive about the dynamics at 
scales larger than that of forcing. 

We wish to emphasize the constancy of fluxes. A constant flux indicates a persistent non-linear 
transfer of energy to smaller scales, i.e. the flow is able to sustain a cascade to arbitrarily small scales, 
regardless of how small the viscous-diffusion parameters are. The term "cascade" necessitates a flux 
constant in wavenumber. It is certainly possible for non-linear interactions to yield a transient transfer 
to smaller scales but one which does not persist (decays to zero) and cannot carry the energy all the way to 
molecular scales. This is sometimes observed, for example, in 2D turbulence simulations and experiments 
(e.g. Figure 1 in [30 ). where we know that a positive downscale flux of energy is only transient and 
cannot be constant (e.g. Figure 1 in |31p. Such distinction between transient and constant fluxes is 



^ The flux in the rS case is noisy at smaU scales because of a highly anisotropic cascade across wavevectors with large 
vertical component k^ while being entirely suppressed across modes with large horizontal component A;^ |15l . The number 
of such anisotropic wavevectors does not increase with K and, therefore, does not provide the additional averaging needed 
to smooth out small-scale noise. 



imperative to modeling efforts. In the former ease, there is no cascade or enhancement of dissipation 
due to turbulence whereas in the latter case, dissipation becomes independent of Reynolds number. 

This issue is especially important when drawing conclusions from limited resolution simulations, true 
of most cases including ours. A constant flux indicates that dissipation should be independent of the 
simulation resolution. One may contend that a constant flux is just a consequence of a steady state and 
having forcing localized to the largest scales. However, it may very well be that the flow reaches steady 
state due to direct viscous dissipation acting on all scales as shown in [5] for potential enstrophy rather 
than a Kolmogorov-like inertial cascade. 
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Figure 2: Plots of total energy flux, (Hk), from runs Rs (top), rS (middle), and RS (bottom). Dotted 
(blue) lines are fluxes taken from an instantaneous time snapshot. Solid (black) lines are averaged over 
time. 

We also compute the potential enstrophy flux. Figure [3] shows that, in a manner similar to that 
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of energy, there is a positive and constant flux of potential enstrophy in all three extreme casea^ The 
plots in Figure [3] are the first measurements of potential enstrophy flux in rotating stratified Boussinesq 
flows and constitute one of the main results in this Letter. They can be regarded as the first empirical 
confirmation of analytical results in [5] which derived an exact law for potential enstrophy flux in physical 
space as a function of scale. Unlike for energy, (fl]^) in Figure k3I for the RS case is not negative over 
K <A. 
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Figure 3: Plots of potential enstrophy flux, (H^), from runs Rs (top), rS (middle), and RS (bottom). 
Dotted (blue) lines are fluxes taken from an instantaneous time snapshot. Solid (black) lines are averaged 
over time. 

Rotating and stratified flows are often said to be 'two-dimensionalized' in some sense, eliciting com- 
parisons with two-dimensional turbulence and often justifying the study of the latter as a simplified 
paradigm for geophysical flows. Here we point out that the existence of a concurrent flux of both energy 



The flux (n^) in the rS case is noisy at small scales for the same reasons as in Figure 2 for {Tlj^). 
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and potential enstrophy in rotating and stratified flow to smaller scales is in itself a marked departure of 
these flows from 2D turbulence. In the latter case, it is known (e.g. [2l|3T]) that the two cascades cannot 
co-exist over the same scale-range since a forward cascade of enstrophy acts as a constraint leading to 
an inverse cascade of energy to larger scales. 

In rotating and stratified flows, the velocity and scalar fields can be decomposed into the sum of a 
vortical component, which accounts for all the potential vorticity in the flow, and into a wave component, 
which has zero potential vorticity (e.g. [271112] ). In the strongly rotating and strongly stratified regime, 
it is known (see [571 HSl IMj) that the vortical component of the flow is governed by quasigeostrophic 
dynamics, which is very similar to 2D turbulence [3]- In this regime, a forward energy cascade of the 
vortical component is suppressed due to the forward cascade of potential enstrophy. We verified (to 
appear in [[15. ) that this is indeed the case in our RS run, where the forward energy cascade (in Figurel2| 
bottom panel) is due to the wave component of the flow in agreement with [27] . However, the situation 
in the remaining two cases, runs Rs and rS, is markedly different from both 2D and quasigeostrophic 
dynamics as will be discussed in J15j . 

7 Conclusions 

The two main results presented in this Letter are (i) energy and potential enstrophy budgets which 
resolve the dynamics simultaneously in space and in scale and (ii) concurrent and persistent cascades of 
energy and potential energy to small scales in three extreme cases of rotating and stratified Boussinesq 
flow simulations. The numerical results on constant fluxes of potential enstrophy constitute the first 
direct empirical evidence in support of analytical results by iBJ. Our findings show a clear departure of 
the flows we study from 2-dimensional turbulence and should be incorporated in any phcnomenological 
treatments of strongly rotating and/or stratified Boussinesq flows. In a longer forthcoming work |15) . 
we shall refine our analysis to study anisotropy of these cascades, their pointwise and scale-locality 
properties, and quantify contributions from vortical and wave components. 
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8 Appendix: Budgets 

For the sake of completion, we write down the complete expression for the transport term in ^ : 

luP luP 

J,(x) = Uj——+U.,T{Ui,Uj)-l^dj^:;-, (16) 






In large-scale potential enstrophy budget ( |ll| , J''(x) is space transport, Z)^(x) is dissipation, and 
^in-i is ^^^ potential enstrophy injected due to external forcing. These are defined as 



-mj 



Df{^) =vN^\Vu:,\'+nf\Vd,W^{v+K)fN{Vd,e).Vu:, (17) 

jf (x) = u, ^-+ q t{u, , q)+ vN q dfj, - nf q d,d,9, (18) 

eZ^ (x) ^ -NqT^, + fq dj^ (19) 



Note that while the first two dissipation terms in (171 are positive definite, the last term can be of either 
sign. 
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